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Abstract. Recent positron flux fraction measurements in cosmic-rays (CR) made by the 
AMS-02 detector confirm and extend the evidence on the existence of a new (yet unknown) 
source of high energy electrons and positrons. We test the gravitino dark matter of bilinear 
R-parity violating supersymmetric models as this electrons/positrons source. Being a long 
lived weak-interacting and spin 3/2 particle, it offers several particularities which makes it 
an attractive dark matter candidate. We compute the electron, positron and 7 -ray fluxes 
produced by each gravitino decay channel as it would be detected at the Earth’s position. 
Combining the flux from the different decay modes we are able to reproduce AMS-02 measure¬ 
ments of the positron fraction, as well as the electron and positron fluxes, with a gravitino 
dark matter mass in the range 1—3 TeV and lifetime of ~ 1.0 — 0.7 x 10 26 s. The high 
statistics measurement of electron and positron fluxes, and the flattening in the behaviour 
of the positron fraction recently found by AMS-02 allow us to determine that the preferred 
gravitino decaying mode by the fit is W ± T Zf , unlike previous analyses. Then we study the 
viability of these scenarios through their implication in 7 -ray observations. For this we use 
the Extragalactic 7 -ray Background recently reported by the Fermi -LAT Collaboration and 
a state-of-the-art model of its known contributors. Based on the 7 -ray analysis we exclude 
the gravitino parameter space which provides an acceptable explanation of the AMS-02 data. 
Therefore, we conclude that the gravitino of bilinear R-parity violating models is ruled out 
as the unique primary source of electrons and positrons needed to explain the rise in the 
positron fraction. 

Keywords: dark matter experiments, cosmic ray experiments, gamma ray experiments, dark 
matter theory 


1 Introduction 


There is vast and compelling evidence of the existence of dark matter (DM), ranging from 
galactic to cosmological scales [1]. Nevertheless, its fundamental nature remains as a mys¬ 
tery [2, 3]. An indirect way to detect DM at Earth could be by measuring its decay or 
annihilation products (i.e. Standard Model (SM) particles), which might be injected in the 
interstellar medium (ISM) by these processes. Thus, the study and characterisation of differ¬ 
ent CR species may shed light on the DM nature. In this respect, the unexpected increase in 
the positron fraction (e + /(e + + e - )) at energies above 10 GeV initially detected by PAMELA 
[4], confirmed by the Fermi-Large Area Telescope (Fermi-LAT) [5], and measured with pre¬ 
cision by the Alpha Magnetic Spectrometer (AMS) [ 6 ] attracted a lot of attention from the 
astro-particles community causing a spate of works trying to explain this anomaly with by¬ 
products of DM annihilation or decay, see e.g. [7-18]. However, the interpretation of these 
leptonic CR measurements in terms of DM is not straightforward as the models of DM able 
to explain the rise in the positron fraction, usually, overpredicts the anti-proton and 7 -ray 
observations, see e.g. [19—24] 1 . Astrophysical hypotheses explaining the positron rise, as for 
instance nearby pulsars or spalls from primary CRs 2 interacting with the ISM, work well re¬ 
producing the data, see e.g. [30, 31]. Nevertheless, these proposals have their own difficulties 
in each case [32], 

The AMS Collaboration has extended the measurement of both, positron fraction and 
positron flux to 500 GeV, as well as the electron spectrum to 700 GeV [33, 34], The main 
highlights of these measurements are: i) the positron fraction rise starts to flatten at energies 
above ~ 200 GeV [33], and ii) precision measurements of electron and positron fluxes allow 
determination of different behaviour of their spectral indexes as the energy increases [34], 
Therefore, in light of the latest AMS-02 release [33, 34] we can better determine the properties 
of the cosmic source responsible for the rise in the positron fraction. 

A minimal model to characterise the electron and positron fluxes, and the anomalous 
behaviour of the positron fraction can be written as [ 6 , 33, 34], 


$ e -(E) = C e -(E/1 GeV) "A- +C S {E/1 GeV)~^e~ E/Es , 

4* e+ (F) = C e +(E/1 GeV) _7e+ + C S [E /1 Ge\)~^ e~ E/Es , (1.1) 

where each equation is composed of an astrophysical background term (modelled with a 
power-law spectrum), and a single source term (modelled with a power-law spectrum with 
an exponential cutoff energy) which is common to the electron and positron flux parameter- 
isations. The astrophysical background term in the electron flux attempt to model primary 
and secondary 3 production, while in the positron flux equation represents only secondary 
production. This minimal model describe in good agreement the positron fraction and the 
combined electron and positron flux [ 6 [. Also, this parametrisation reveals and characterises 

x Note that decaying DM particles could affect the physics of dense stars implying strong constraints on the 
DM properties [25]. These constraints are avoided by the gravitino because its elastic scattering cross section 
on nuclei is proportional to 1/Mpi [26], therefore there is no chance that they can be trapped in the core of 
dense stars. In particular, antideuterons produced by gravitino decays were recently studied in [27-29]. 

2 The most common primary CR particle is the ubiquitous proton or hydrogen nucleus. 95% of all CRs 
are protons, 4% are helium nuclei, and the 1% balance is made up of nuclei from other stellar-synthesized 
elements up to iron. Their principal Galactic source are: supernova remnants, neutrons stars, black holes and 
possibly DM 

j It is due to spallation in the ISM owing to the impacts of primary cosmic rays off ISM nuclei. 
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the spectrum and intensity of a common source of primary electrons and positrons required 
to fit the high energy end of the spectra. 

In this work we assume the single common source term in eq. 1.1 consisting of decaying 
DM particles, which preserves CP symmetry, and thus contribute an equal amount to the 
electron and positron flux. In particular, we study a decaying gravitino DM scenario, which 
is possible in the context of R-parity violating (RpV) supersymmetric (SUSY) models, see 
e.g. [21, 35-43] 4 5 . We assume that the relic density produced by our gravitino DM, is consistent 
with the observations, as already has been shown in several places [44-48]. It is worth noting 
that these works assume that the gravitino is sufficiently long-lived in order to maintain its 
comoving density from the period of decoupling until the present. In fact, since RpV terms 
are in general quite small and the gravitino interactions are suppressed by the Planck mass, 
the latter condition is naturally obtained [36]. 

In SUSY with RpV, the gravitino decay channels are fixed (see for instance [35-37, 49, 
50]). In this work we just consider the primary two body decay channels which could pro¬ 
duce electrons and positrons in the final state. In models with trilinear RpV two body decay 
channels are possible at one-loop level, but for heavy gravitinos the behaviour of the total 
decay rate is dominated by three-body decay channels at tree level [42, 43]. Therefore, we 
implicitly assume negligible R-parity violation due to trilinear couplings, so our theoretical 
framework are the bilinear RpV models. For several values of the gravitino mass, final grav¬ 
itino decay products are computed using Pythia 8.185 [51] and propagated from regions with 
non-zero DM densities towards the Earth. The obtained spectrum produced by each of the 
different decay channels are combined in order to fit the most recent AMS-02 data. From 
the results of this fitting procedure, we can extract the best values for the corresponding 
branching ratios (BR). We then compute the corresponding diffuse 7 -ray emission produced 
by the obtained gravitino parameters at high Galactic latitudes, which turns out to signifi¬ 
cantly exceeds the flux of the extragalactic 7 -ray background (EGB) recently measured using 
the Fermi-LAT [52] . Finally, we compute the residual upper limit flux of the EGB removing 
its known contributors, in particular the integrated emission from blazars as modelled in [53], 
and then compare it to the prompt gravitino-induced 7 -ray flux to obtain conservative limits 
on the gravitino lifetime that can be directly compared to the lifetime required to account for 
the AMS-02 leptonic measurements. 

This work extends and update previous studies in several respects. As we try to be 
as model independent as possible, we are not fixing the branching ratios for the different 
gravitino decay channels, neither the flavour structure in our theory, but rather, this is the 
outcome of the fits to the AMS-02 data. Indeed, we present for the first time a fit of the 
branching ratios for the different gravitino decay channels to data, finding that the best fit 
parameters are compatible with the theoretical expectation for gravitinos with masses in the 
range between 750 and 3000 GeV, in proportion 2:1:1 [22, 26, 50]. Further¬ 

more, we are using the most up to date results from AMS-02 and Ferrm-LAT collaborations, 
which considerably extends previous measurements. It allows a more precise characterisation 
of the source responsible for the positron fraction excess and reduce the window for exotic 

4 Unlike [36, 38-41] we focus on gravitino DM with mass above 1 TeV motivated by the flattening of the 
positron fraction spectrum at ~ 200 GeV. The study of the consistency between this gravitino mass scale and 
other observables, as neutrino physics, is beyond the scope of this work. 

5 The EGB comprises detected point sources, all extragalactic emissions too faint or too diffuse to be 
resolved, and any residual Galactic foregrounds that are approximately isotropic. Therefore, the EGB could 
contains emission from gravitino DM decays of extragalactic and Galactic origin, it has been used to set 
stringent constrains on DM lifetime, see [23] and references therein. 
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contributors to the 7 -ray sky. In particular, the flattening in the behaviour of the positron 
fraction at about 200 GeV and the high statistics measurement of electron and positron fluxes 
allow us to determine the mass scale of the gravitino DM (~ 1 TeV ) and that the gravitino 
decaying mode preferred by AMS-02 data is W ± r =F . The last differs from previous analyses 
where the favoured channel is IT ± /i T [ 20 , 21 ], we discuss further on this difference in section 
3.2. It is worth noting that we only use in the fit the positron fraction and positron flux 
measurements in order to increases the reduced \' 2 as the number of degrees of freedom is re¬ 
duced. We do this to emulate the effect of assuming a covariance matrix correlating the three 
datasets, which is beyond the scope of this work. We ruled out the hypothesis of gravitino 
DM in bilinear RpV models as the sole explanation for the rise in the positron fraction using 
the latest determination of the EGB [52], 

The paper is organised as follows: In section 2 we introduce the primary gravitino decay 
channels which are available in SUSY with bilinear RpV models. In section 3 we consider the 
propagation of electrons and positrons produced by the gravitino decays through the ISM in 
three different propagation model scenarios. Afterwards we fit the most recent leptonic data 
published by the AMS Collaboration using a model that consider a power-law background 
plus a source given by the gravitino decays 6 . In section 4 we use the results obtained from the 
fit to AMS-02 data in order to determine the 7 -ray signal from Galactic gravitino DM decay 
at high latitudes and we confront these results with the new Fermz-LAT EGB measurement 
and latest modelling of integrated blazar emission. In section 5 we use the EGB measurement 
to constrain the gravitino DM models that provide an acceptable explanation to the positron 
rise. Finally, in section 6 we conclude. 

2 Gravitino decay modes 

Our decaying DM candidate consists of a heavy gravitino that decay through RpV terms. As 
already pointed out, we will study gravitino LSP. We consider m 3 / 2 ~ 1 TeV motivated by 
the most recent positron fraction measurement presented by the AMS Collaboration. In this 
regime the main dacay modes of the gravitino are the following [37, 50], 

T3/2 -> W±lT, Zv u (2.1) 

where i is the lepton family index'. It has been shown in [22] that in the limit of gravitino 
masses above 1 TeV, the decay widths of the channels listed in 2.1, reach their asymptotic 

6 As already mentioned, we only use positron flux and positron fraction measurements in the fit, then we 
use the obtained gravitino and background parameters to check visually a reasonable agreement of the electron 
flux with the measurement. 

'The Gravitino decay to 717 is strongly suppressed in the limit where m 3 / 2 is heavier than the weak 
gauge bosons and will not be further considered in this work, also we do not consider decays of the gravitino 
involving the extra Higgs bosons predicted in SUSY theories, because they are generally heavier than the 
LSP [22, 26, 54], 
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Where Aj are the effective bilinear RpV terms [55], fi is the higgsino mass and v the vacuum 
expectation value. The corresponding branching ratios become independent of the size of 
R-parity violation, while their relative sizes converge to the proportion 2:1:1, for , Z zv ? ;. 

and Hui, respectively. From now on we generically refer to these final states as AZj. Since 
we are indeed interested in the heavy gravitino regime, where the hierarchy of the branching 
ratios for the different decay channels holds, we assume this as our benchmark scenario. 

We use Pythia 8.185 [51] to compute the number of final state electrons, positrons and 
photons, as well as their 4-momenta, produced by each single gravitino decay mode. In fig. 1 
we show the resulting energy spectra. 


m 3 / 2 = 2 TeV, positrons m 3/2 = 2 TeV, photons 




Figure 1 . The positron (a) and photon (b) energy spectrum produced by the different gravitino 
decay modes as a function of energy. The electron distribution is equivalent to the positron case. 


3 AMS-02 positron excess and gravitino dark matter 
3.1 Cosmic ray propagation 

Once a gravitino decay in the halo of the Milky Way (MW), some of the resulting charged 
particles will propagate through the ISM towards the Earth, where they might be detected. 
Theoretically, the propagation of electrons and positrons in the MW can be treated as a 
diffusive process [56]. If we assume that this is a steady state process, we can estimate the 
electron and positron density per unit energy at position r (in Galactic coordinates) i.e. 
f(E,r), by solving the following diffusion equation, 
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(3.1) 


o 

- V • ( JC(E , r) V/) - — (b(E, r)f) = Q(E, r), 

where JC(E,r) = ICo(E/GeV) s is the diffusion coefficient that account for the transport 
of electrons and positrons through the Galactic magnetic fields, b(E, r ) is the energy loss 
coefficient due to the different types of interactions of electrons and positrons with the ISM, 
like synchrotron radiation, bremsstrahlung and inverse Compton scattering (ICS), and Q is 
the source term which model how electrons and positrons are injected in the ISM, in our case, 
this is due to gravitino decays, and is given for each A/* channel as: 


Q{E,r) 


P(r) \ 1 dN \i^ e ± 

"*3/2/ r 3/2 dE 


(3.2) 


Replacing this expression for Q(E,r) in eq. 3.1, we can determine the flux of elec- 

o /o 

trons and positrons as a function of the energy E and position r, i.e. d$>'^ ^ e ±/dE ( E,r) = 
v e ±f(E,r)/4iT, using the following equation [57] 8 : 
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dE s 
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3/2 


A/j— 
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(E s ) I(E,E s ,r), (3.3) 


where v e ± is the velocity of the electrons and positrons, E s is their energy at the source and 
I(E, E s , r) are the generalized halo fimctions (GHF), which are essentially the Green functions 
connecting a source with a fixed energy E s to a given energy E (see [58] for more details of 
the GHF). In this work we use the Navarro, Frenk and White (NFW) [59] density profile 


Pnfw(^) = Ps — 
r 



(3.4) 


where we adopt r s = 24.42 kpc and p s = 0.184GeVcm~ 3 following [57], where at Earth 
position 7 *q = 8.33 kpc [60] the local DM density is assumed to be p@ = 0.3 GeV cm" 3 . 
The diffusion model is defined by /Co, 6 and the half-thickness L of the diffusion zone 9 . We 
use the set of propagation models called: MIN, MED, MAX described in [58]. The MIN 
(MAX) model minimizes (maximizes) the p flux at Earth and are in agreement with the 
measurements of the boron-to-carbon (B/C) ratio at the Earth’s position [58]. This models 
contains parameters for the propagation model that achieve an acceptable configuration for 
CRs propagation in the MW. The use of these models are aimed at estimating the impact 
of the uncertainties inherent to electrons and positrons propagating in the ISM. Note that 
MIN, MED, and MAX scenarios do not test the complete range of systematics related to the 
modelling of CRs propagation [58]. Furthermore, the MIN scenario is strongly disfavoured 
by CR and gamma-ray data [61, 62], 

In panel a) of fig. 2 we show the propagated electron spectra at the Earth’s position for 
different gravitino masses using MED propagation model. We show the propagated spectra 
of electrons from gravitino decays with masses of 1, 2 and 3 TeV. We use this set of curves in 
the following section as the inputs to the fit of the AMS-02 data. In panel b) of figure 2 we 


8 In [58] the eq. (3.1) is resolved for the MW configuration. 

9 It is define as a solid flat cylinder with height 2 L in the z direction and radius R = 20kpc in the radial 
direction. 
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Figure 2. a) Propagated flux at Earth using the MED scenario for a gravitino mass of 1 TeV, 2 TeV 
and 3 TeV with fixed 73/2 = 10 26 s, different colours indicates the different channels of the gravitino 
decay, and we are assuming only one decay mode at a time, with a BR\i t equal to one. b) Propagated 
flux at Earth using the MIN, MED, MAX scenarios for a gravitino mass of 2 TeV. We note that MIN 
scenario is strongly disfavoured by recent CR and gamma-ray data [61, 62]. It is worth noting that 
MIN and MAX scenarios do not bracket all the uncertainty due to propagation of CRs in the ISM. 


present the effect on the gravitino-induced electron flux at Earth’s position due to the three 
different propagation models MIN, MED, and MAX. We see in this figure that the effect of 
different propagation scenarios is less than one order of magnitude, but in any case as we will 
see later the results of the fit are barely affected by the choose of the propagation model. 

3.2 Determination of gravitino parameters from AMS-02 data 

In this section we consider a decaying gravitino scenario as describing the anomalous be¬ 
haviours of the electron and positron fluxes, and the positron fraction, all of them recently 
measured with high precision by the AMS-02 experiment [33, 34], Since the three mea¬ 
surements are not fully independent among them, correlations should be included in the fit, 
for instance by using a complete covariance matrix. Unfortunately, the computation of this 
matrix is beyond the scope of this work. Instead of a full statistical analysis assuming a 
covariance matrix we make an approximation in the fit computation. In practice, we just 
take into account two of the three datasets plus a diagonal covariance matrix, which can be 
computed using the standard deviations released by the AMS Collaboration. In this way we 
reduce the artificial over constraining of the free parameters that we would obtain from the 
fit using the three data sets disregarding correlations among them. 

From the three possible combinations of datasets we consider the positron flux and 
positron fraction, excluding the electron flux. By using these datasets it is possible to deter¬ 
mine unambiguously both, the best fit values of the free parameters and their error intervals. 
Similar results, both qualitatively and quantitatively, are obtained at using the electron flux 
and positron fraction datasets. For the remaining case, electron and positron fluxes, it turns 
out that the normalisation and spectral index of positron flux, C e + and 7 e + , cannot be deter¬ 
mined. It is so because the high energy part of the measured positron flux can be well fitted 
alone by the gravitino-induced positron emission, therefore we just need a mild contribution 
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of the positron flux at low energies that can be arranged by using unconstrained combinations 
of C e + and 7 e+ . 

Motivated by eq. 1.1 we model the electron and positron flux by considering the sum 
of a non-symmetric power law background plus the contribution of a symmetric source, in 
our case, produced by gravitino decays in the the Galactic halo, which is assumed to make 
the whole DM density in the MW. The latter assumption can be relaxed when we consider 
that a different value of the average density can be absorbed by a redefinition of the gravitino 
lifetime. The only constraint would come from the metastability condition, that roughly 
speaking requires r 3 / 2 > 10 17 s. The electron and positron fluxes are given by, 

®e-(E) = C e -E~^- +<f> 3 e /2 (E,m 3/2 ,T 3/2 ,BR Xh ) : 

$ e+ {E) = C e+ E-^ +$H\E,m 3/2 ,T 3/2 ,BR xli ). (3.5) 

The asymmetric contribution to the eE flux is determined by the corresponding normali- 

o /o 

sation factors C e T and spectral indexes 7 eT . The term $ ^ ( E , m 3 / 2 , T3/2, BR\i t ) corresponds 
to the symmetric contribution to the flux arising from the decay of gravitinos. As already 
mentioned in section 3.1, the gravitino contribution depends on its mass m 3 / 2 , lifetime r 3 / 2 
and decay branching ratios BR\ Notice that one of the gravitino BRs is not independent 
because we have assumed that the sum of them must be equal to unity. Apart from the latter 
restriction, during the fit we consider these variables as free parameters. 

Since the computation of the gravitino contribution to the eE flux is quite sensitive to the 
value of its mass, for simplicity we have considered just three specific cases given by = 
1, 2 and 3 TeV. We will see later that these cases are enough to extract general conclusions 
from the fit to the data. For each of these masses we compute a set of spectra corresponding 
to the different decay channels with fixed lifetime r.^ 2 = 10 26 s after propagation through the 
ISM. These spectra are showed in fig. 2. The flux obtained for this specific lifetime can be 
linearly scaled in order to obtain the flux for any other value of r 3 / 2 . Thus, we model the 
gravitino contribution in the fit computation with the following expression, 

$%\E,m* 3/2 ,T 3/2 ,BR xli ) = 3 /2 >Tb* /2 )- (3.6) 

3/2 A U 

Finally, in order to determine the best fit values of the free parameters, for each of 
the rny 2 considered we minimise a %-squared function, x 2 (C e +, 7 e +, C e - , 7 e -, r 3 / 2 , BRx^), by 
considering the eq. 3.5 and the data points and systematic errors from positron flux and 
positron fraction [33, 63]. Indeed, we sum two x-squared functions constructed from each 
data set. We use 48 points from positron flux (E > 10 GeV) plus 65 points from positron 
fraction (E > 1 GeV). 

In table 1 we report the values of the best fit parameters using the reference MED 
propagation model. Also, we compute the error intervals at la confidence level for the free 
parameters following [64], It is worth noting that the lifetimes obtained in the fits are close 
to the nominal value Ty 2 = 10 26 s in the three scenarios, while the dominant BRs are given 
by the channels T 3 / 2 —> W ± T :f , T 3 / 2 —> Zv and T 3/ / 2 —> VF ± e =F (^) =F in descending order. 

Including the theoretical prediction for the BR proportions 2:1:1 (see section 2) in the 
fit we obtain reasonable p-values (p>0.05), but the reduced x-squared is quite large in com- 
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parison to the minimum obtained with free BRs 10 . The resulting gravitino lifetimes when the 
theoretical prediction on the BRs is included are lower than in the case with free BRs, they 
are in the range 5.9 — 6.5 x 10 25 s, these lifetimes are excluded by the limits from antiproton 
measurements derived in [22], The antiproton limits in [22] cannot be directly applied to our 
results with free BRs, which provide better description of the AMS-02 leptonic data, therefore 
in next sections we use gamma-ray measurements to probe those gravitino scenarios. 

For the fit with free BRs we find that the main gravitino decaying channel is W ± r T . 
This result differs from conclusions in [20] and [21] where the gravitino decay mainly into 
the second lepton generation. In the following we go further on the differences between those 
analysis and this work that make the channel the preferred by the fit to explain the 

AMS-02 leptonic data. In [20] the PAMELA positron fraction [4] and Fermi -LAT electron 
plus positron flux [65] are interpreted in terms of decaying DM. In particular one important 
difference between the data used in [20] and this work is that the positron fraction measured 
by PAMELA is harder with respect to the fraction measured by AMS-02, see figure 3. It 
favours lighter families, as can be seen in figure 2; i.e. for lighter families, both, electron 
and positron fluxes become harder. Also, it is important to note that [20] only studies pure 
DM decaying channels unlike our approach where we combine all of them. The other work 
[21] concludes that the W is the favoured gravitino decaying channel by comparing the 
positron fraction yield by a benchmark point using two different pure decay channels, W^eF 
and VF ± /T F , to the AMS-02 data. In [21] there is no fit to AMS-02 data, and the electron 
and positron background fluxes are determined using the electron flux measured by PAMELA 
[66]. Therefore, the inclusion of high statistic new data and the reduction of error bars provide 
us new insight in the determination of the source responsible for the positron fraction rise, 
making plausible our new result. 

In fig. 3 the positron fraction obtained from the best fit values is compared to the 
PAMELA and AMS-02 data, while a similar comparison for the electron and positron flux 
are shown in fig. 4. Although we do not consider the electron flux measurement in the fit we 
use the parameters obtained in the fit to compute the electron flux. Notice that a gravitino 
with 7713/2 = 1 TeV fails to adjust the highest data point of the electron flux spectrum reported 
by AMS-02 (see left panel of fig. 4), which therefore disfavours scenarios with 7713/2 < ITeV. 
From the behaviour of the lines corresponding to 7713/2 = 2 TeV and 3 TeV we can see that we 
are at the edge of the measured spectra, and that any higher mass value will over predict the 
data at high energies. 

For completeness we repeated the fit using MIN and MAX propagation models. We find 
for the three gravitino masses under study that the best fit parameters are compatible with 
results reported in table 1. 

4 The Extragalactic Gamma-ray Background measured by Fermi -LAT and 
gravitino dark matter 

The 7 -ray sky has been observed by the Fermi -LAT Telescope with unprecedented detail. 
Most of the 7 rays detected come from: i) point-like or small extended sources and ii) a 
strong diffuse emission correlated with Galactic structures [67]. In addition, a tenuous diffuse 
component has been detected, the Isotropic 7 -ray Background (IGRB) [ 68 ]. The origin of the 
IGRB can be sources that remains below the detection threshold of the Fermi-LAT, among 

10 The Ax 2 with respect to the minimum turns out to be larger than 11.78, which determines the la region 
around the global minimum. 
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Figure 3. The positron fraction measured by AMS-02 and PAMELA (blue and cyan points respec¬ 
tively), and fits using the gravitino decay products + background for three different values of the 
gravitino mass (in red), 1 (continuum line), 2 (dot-dashed) and 3 TeV (dashed). The background and 
gravitino contributions are also shown separately in grey and pale brown respectively. 




(a) (b) 

Figure 4. The electron (a) and positron (b) flux measured by AMS-02 and PAMELA (blue and cyan 
points respectively), and fits using the gravitino decay products + background for three different values 
of the gravitino mass (in red), 1 TeV (continuum line), 2 TeV (dot-dashed) and 3 TeV (dashed). The 
background and gravitino contributions are also shown separately in grey and pale brown respectively. 

others, DM decay/annihilation [69, 70] can produce a sizeable contribution 11 . The observed 

“The main EGB contributors are blazars, star-forming galaxies, diffuse processes as intergalactic shocks [71- 
73], interactions of ultra high energy CRs with the Extragalactic Background Light (EBL) [74], and CR 
interactions in small solar system bodies [75] 
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41.49 

47.26 

35.88 

36.22 

36.4 

7e+ 

3.82 

3.89 

3.8 

3.86 

3.78 

3.84 

3.78 

3.74 

3.73 

73/2 [10 26 s] 

0.99 

1.04 

0.79 

0.89 

0.68 

0.79 

0.78 

0.66 

0.59 

BR w ± e ^ 

0.09 

0.12 

0.05 

0.11 

0.0 

0.13 

0.05 

0.0 

0.0 

BR\v±fiT 

0.0 

0.08 

0.0 

0.19 

0.07 

0.24 

0.0 

0.0 

0.0 

BR w ± t 

0.91 

0.93 

0.82 

1.0 

0.68 

1.0 

0.61 

0.44 

0.32 

BRzv 

0.0 

0.29 

0.12 

0.40 

0.25 

0.47 

0.0 

0.0 

0.0 

BR*Hu 

0.0 

0.23 

0.0 

0.11 

0.0 

0.13 

0.0 

0.0 

0.0 


Table 1. Results of AMS fit in a decaying gravitino scenario using MED propagation model. During 
the fit we have considered as the dependent branching ratio, however in order to compute their 

confidence intervals we consider BRz v as the dependent one. We use N = 104 as the effective number 
of degree of freedom. The corresponding p-values for each scenario are equal to 0.99. 


IGRB depends on the point source detection threshold of the instrument. A physical quantity 
is the total EGB, defined as the combination of resolved sources and the IGRB. The Fermi- 
LAT Collaboration has a new determination of the EGB from 50 months of data that expands 
from 100 MeV to 820 GeV, in fig. 5 we present this measurement (cyan points with error 
bars) [52] 12 . In the extraction of the EGB a relevant source of systematic uncertainty is the 
modelling of the Galactic foreground (FG) emission 13 , the yellow band in fig. 5 represents 
this uncertainty. 

In fig. 5 in addition to the Fermi -LAT EGB (cyan points with error bars [52]) we present 
the 7 -ray emission due to known EGB contributors, star-forming galaxies (red band [76]) 
and radio galaxies (blue band [77]) as well as the integrated emission of blazars with EBL 
absorption as recently modelled in [53] (green band). The grey band represents the sum of 
all these components, it accounts for the observed amplitude and spectral shape of the EGB. 

12 This is the EGB measurement using FG model A. 

13 In order to test this systematic, the fitting procedure to obtain the EGB was applied to many different 
Galactic FG models. In particular, reference FG model A is derived from the class of models presented in [68]. 
Models B and C have more degrees of freedom in electron and diffusion coefficient, for details on FG models 
A, B and C we refer to [52]. 
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From this figure we can see that there is little room for other contributors, such as DM 14 . 
We compute the 95% CL upper limit flux subtracting the average emission from non-exotic 
contributors (black line) to the EGB (cyan points) and adding 1.64 times the sum of data 
and model uncertainties in quadratures 15 . 


4.1 7 -ray flux at Earth from gravitino decays 

EGB contributions from gravitino DM decays can has 3 different sources: i) the smooth 
Galactic halo, ii) sub haloes hosted by the Galactic halo and iii) extragalactic structures. The 
signal from extragalactic gravitino decay is expected to be isotropic, although it is independent 
on the amount of DM clustering at each given redshift [79-81], it is subject of attenuation 
which increase significantly with the y-ray energy [82], The emission from the smooth Galactic 
halo provides a lower limit on the y-ray total gravitino contribution to the EGB, we only use 
this emission in the analysis. 

We calculate the differential flux of y rays from gravitino decays in the Galactic halo 
by integrating the DM distribution around us along the line of sight. We consider gravitino 
decays producing y rays in the final state. In this case the flux reads: 


d4>J) al ° 

dEdn 


1 

47 tt 3 /2 rn 3 / 2 


1 d./V lotal 
AD dE 


/ An 


cos bdbd£ / ds /9haio( r (s, 5, £)), 

Jo 


(4.1) 


where b and £ denote the Galactic latitude and longitude, respectively, and s denotes the 
distance from the Solar System. Furthermore, AD is the region of interest (ROI). The radius 
r in the DM halo density profile of the Milky Way, phalo is expressed in terms of these Galactic 
coordinates 

r(s, b, £) = 's 2 + 7q — 2 sr Q cos b cos £ . (4.2) 


For gravitino masses above the Z boson mass, the total number of photons produced in 
gravitino decays can be expressed as 


dlV lotal 

dE 


\h 


dN 


3/2 


Mi 


dE 


(4.3) 


o /o 

where dN^ ^/dE is the photon energy spectrum produced by different gravitino decay 
channels computed in section 2, see panel b) of fig. 1. 

In the next section we contrast the y-ray flux predicted by the gravitino DM that can fit 
the AMS-02 measurements with the latest EGB measurement with the Fermi -LAT Telescope. 


5 Contrasting the rise in positron fraction with the EGB through gravitino 
dark matter 

In fig. 5 we show the contributions to the EGB from the three gravitino decay scenarios 
presented in table 1, which provide good fits to the AMS-02 data (gravitino signal red line). 
We also plot the gravitino-induced gamma-ray emission using the results when we include the 

14 Note that Galactic templates of DM decay can be degenerate with ICS templates of the Galactic FG 
emission model as pointed out in [78]. Nevertheless, we only contrast the gravitino y-ray emission to the EGB 
in an energy range (> 500 GeV) where the ICS component is not very important, as shown in [78]. 

15 We assume Gaussian errors, therefore 95% of the area of a Gaussian distribution is within 1.64 standard 
deviations of the mean. 
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Figure 5. We contrast the predicted 7 -ray flux in the three scenarios obtained from the fit to AMS-02 
data to the EGB as measured by Fermi- LAT, 7713/2 = 1 TeV, r 3 / 2 = 0.99 x 10 26 s; m 3 / 2 = 2 TeV, 
T3/2 = 0.78 x 10 26 s; and m 3 / 2 = 1 TeV, r 3 / 2 = 0.67 x 10 26 s; solid, dot-dashed, dashed dark-red 
lines, respectively. The gravitino DM-induced 7 -ray emission clearly exceeds the EGB limits. The 
EGB as measured by Fermi -LAT using FG model A (cyan points) and its uncertainty due to Galactic 
foreground modelling (yellow band). Brown points represent 95% CL upper limit flux after subtracting 
from the EGB the average contribution (grey band) from star-forming galaxies (red band [76]) and 
radio galaxies (blue band [77]) as well as the integrated emission of blazars with EBL absorption as 
recently modelled in [53] (green band). 


2:1:1 proportion on the BR in the fit (gravitino signal blue line). They all clearly overshoot 
the EGB measured by Fermi- LAT. We already knew that the scenarios with the restriction 
2:1:1 on the BRs were excluded because the antiproton limits derived in [22], On the one 
hand, the results inferred from Fig. 5 show that the EGB limits can be as sensitive as the 
antiproton measurements in order to test this gravitino scenario. On the other hand, the 
antiproton limits in [22] cannot be directly applied to the gravitino scenarios with free BRs. 

As already mentioned, we can see from Fig. 5 that the points giving the best fits to 
the AMS-02 data clearly overshot the EGB measurement. Now, in order to be more general, 
we extend this analysis considering the scenarios that allow a relatively good fit of AMS- 
02. In this way we check whether we can have a gravitino-induced gamma ray emission 
consistent with the EGB limits at the cost of a poorer fit to the positron flux and positron 
fraction measured by AMS-02. In practice, we test points of the parameter space that are not 
excluded at 95% CL by the data, i.e. with a p-value > 0.05. We perform the following analysis 
to check this out. For fixed gravitino mass and lifetime, we call this lifetime T[ e +f rac e +] 16 , we 
marginalise over the BRs and background parameters, with and without their theoretical 
prediction 2:1:1, with respect to the positron flux and positron fraction measurements. We 

16 to point out that comes from an analysis using positron flux and positron fraction. 
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Figure 6. Comparison between the gravitino lifetimes required to provide acceptable fit (p-value 
> 0.05, color code) to positron flux and positron fraction T[ e +f rac ,e+] and the lifetimes required to 
do not overshoot the EGB residual emission, tegB) for different gravitino masses. Points with ratio 
T [e+frac,e+]/ T EGB lower than 1 are excluded since a gravitino with 'T[ e +f rac e +] would produce brighter 
gamma-ray emission than the allowed maximum produced by similar gravitino (same mass and BR 
distribution) characterised by tegb- Squares represent the points with p-value>0.05 when the restric¬ 
tion 2:1:1 in the BR proportions was considered in the fit, the scenarios with free BRs are represented 
by circles. All the points reside below the limit T[ e +f rac e +]/rEGB=l) are excluded. The smaller plot 
in the right-upper corner is a zoom out of the bigger plot which points out the size of the exclusion. 


keep the BR distribution corresponding to the largest p-value, then we repeat the procedure 
for different masses and rr e + f rac ,e+]) i n particular we explore 5 gravitino masses, 1.0 TeV, 
1.4 TeV, 2.0 TeV, 2.4 TeV, and 3.0 TeV. Next, we use the BR distributions obtained for 
each of the masses and T[ e +f rac e +i to compute the minimum gravitino lifetime allowed by the 
EGB residual component (brown points in fig. 5), tegb- In this way for each single point 
in the parameter space [mass, Tr e +f rac?e +], BRs] with p-value > 0.05, we obtain the minimum 
tegb which determine if the point [mass, T[ e + f raCje +i, BRs] is excluded or not, as the ratio 
t[ e+frac,e+] /tEGB is lower or larger than unity, respectively. Results are presented in figure 6, 
where we contrast each of the gravitino DM models that provide an acceptable explanation 
to the positron rise with the limits on gravitino lifetime from the EGB residual component. 
All the points reside below the limit T[ e +f rac e +i/TEGB —1) implying that all the gravitino DM 
scenarios providing relatively good fit to AMS-02 leptonic data are excluded, disregarding 
their BR distribution or mass. 

Our results are in agreement with the limits for gravitino DM presented in [24]. It is 
worth noting that the limits presented in fig. 6 are obtained contrasting the gravitino-induced 
7 -ray emission from the Galactic smooth halo, excluding the Galactic plane region, to the 
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EGB, for a variety mixture of gravitino decaying channels as required to reproduce AMS-02 
leptonic data. This is a different approach to what was used in [24], where the gravitino- 
induced gamma-ray emission considered is originated in extragalactic structures and only a 
single gravitino decay channel is consider W±. As already mentioned, only considering 
the Galactic smooth halo contribution provides a minimal contribution to the EGB emission 
from gravitino decay, but it is enough for the purposes of this work. 

6 Conclusions 

Recent AMS-02 results provide new insights on the anomalous rise in positron fraction, pre¬ 
senting the first evidence of a flattening spectrum at higher energies [33]. In addition, the 
measurement of electron and positron fluxes allow a better characterisation of the phenomenon 
behind the positron fraction rise [34], We model positron and electron fluxes as decreasing- 
with-energy power laws plus an extra source injecting equally electrons and positrons in the 
ISM, we test the gravitino of RpV SUSY models as the source term (e.g., [36-41]). We com¬ 
pute the electron, positron and 7 -ray spectra produced by gravitino DM decays as would be 
detected in the Earth. We have shown that the gravitino DM with mass in the range 1 — 3 TeV 
and lifetime of ~ 1 —0.7 X 10 26 s, mainly decaying to can reproduce the energy spectra 

of the positron fraction and the electron and positron fluxes at the same time. In order to test 
the viability of these scenarios we study the corresponding 7 -ray emission and we find that 
these scenarios clearly exceed the observed EGB. Then, we use the most recent Fermi-L AT 
EGB measurement [52] and a model of its known contributors, star-forming galaxies, radio 
galaxies, and blazars 1 ' to compare with the emission from every single gravitino DM able to 
provide an acceptable explanation of the AMS-02 leptonic measurements, characterised by 
their p-value larger than 0.05. We find all the points overshoot the EGB limits by about one 
order of magnitud, this excludes the possibility of gravitino DM in bilinear RpV models as 
unique source of primary positrons 18 . Thus, if those gravitinos made up the whole DM of the 
Universe, they can only contribute modestly to the highest energetic electron and positron 
fluxes detected at the Earth. 

Acknowledgments 

The authors are thankful to Andrea Albert, Borut Bajc, Marco Cirelli, Michael Grefe, Luis Labarga, 
Carlos Munoz, Paolo Panci, Frank Steffen, and Gabriela Zaharijas for useful comments, and Marco 
Ajello for providing the EGB model contributions of figure 6. Also we thank the anonymous ref¬ 
eree of PDMU for recommendations in the fitting procedure. This work was supported by Conicyt 
Anillo grant ACT1102. GAGV thanks for the support of the Spanish MINECO’s Consolider-Ingenio 
2010 Programme under grant MultiDark CSD2009-00064 also the partial support by MINECO under 
grant FPA2012-34694. BP also thanks for the support of the State of Sao Paulo Research Foundation 
(FAPESP). The work of NV was supported by CONICYT FONDECYT/POSTDOCTORADO/3140559 


11 The integrated emission of blazars include the EBL absorption we use the model in [53] 
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